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(57) A system and methods for indicating formation resistivity anisotropy and evaluating the vertical 
resistivity component by combining array induction measurements and array laterotog measurements. 
Formation measurements are obtained with an Induction tool and a iaterolog tool. Formation parameter data 
are derived from the induction tool measurements. The Iaterolog tool response is modeled using the derived 
formation parameter data. The formation resistivity anisotropy and/or the vertical resistivity component are 
determined from a comparison of the modeled Iaterolog tool response against the apparent measurements 
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METHOD AND SYSTEM FOR INDICATING ANISOTROPIC RESISTIVITY IN 

AN EARTH FORMATION 

1. BACKGROUND OF THE INVENTION 

1.1 Field of the Invention 

The present invention is generally directed to the analysis of underground earth 
formations, and, more particularly, to the determination of formation resistivity properties 
and/or profiles. 

IJl Description of Related Art 

Knowledge of the resistivity properties of underground earth formations is 
fundamental in the evaluation and characterization of potential and existing hydrocarbon- 
bearing reservoirs. The terms "logging" and "log" are respectively used herein to 
designate any sequential measurements and the recording of such measurements as a 
function of depth and/or time of variations in a given characteristic of die formations 
around a borehole. 

Formation resistivity anisotropy is a characteristic of reservoirs that can complicate 
formation evaluation. Many reservoir rocks exhibit resistivity anisotropy, especially when 
saturated with oil. There are several mechanisms, which can produce tiiis anisotropy, 
among which are very thin sand-shale laminations, depositional changes in clean 
sandstone, and wind-distributed sands (aeolian formations). See Rubin, D. M., Cross 
bedding, bedforms, and paleocurrents, Scx:iETY OF ECONOMIC PALEONTOLOGISTS AND 
Mineralogists, Concepts in Sedimentology and Paleontology, 1; Klein et al.. The 
petrophysks of electrically anisotropic reservoirs. Transactions of the SPWLA Thirty- 
Sdcih Annual Logging Symposium, Paris, France, June 26-29, 1995, paper HH. 
Resistivity measurements from witiiin a borehole drilled through these formations have 
been made by a number of techniques employing tht use of well-known instmments, 
including induction and laterolog (electrode) type devices. 

In induction logging, a transmitter coil mounted on a well tool is lowered into a 
borehole and energized by an alternating current Indications are then obtained of the 
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influence of surrounding formations on the electromagnetic field established by the coil. 
Usually such indications are obtained by observing the voltage induced in a receiver coil 
also mounted on the well tool in coaxial relationship with the titmsmitter coil and axially 
spaced apart therefrom. Conventional induction well tools include at least one transmitter 
and a receiver coil array. Such an induction tool is illustrated and described in U.S. Pat. 
No. 5,157,605, assigned to the present assignee. 

The laterolog technique comprises making a measurement with a well tool having a 
current emitting electrode and voltage electrodes one either side thereof which are 
operated to force current into the formation. The current is focused into the formations 
and returns locally. i.e., to the tool body. The formation resistivity is determined by 
measuring the current flow. Conventional laterolog tools are illustrated and described in 
U.S. Pats. Nos. 5.754,050, 5,852,363, 6,023.168, and Documents EP-478409, EP-544583, 
all assigned to the present assignee. Docmnent EP-544583 describes a laterolog tool 
including botii a "deep" investigation mode and a "shallow" (LLS) investigation mode. 

The sensitivity of a resistivity logging tool to resistivity anisotropy depends on the 
physics of tiie measurement and tiie relative dip angle of the well bore of the formation. 
For example, 2MHz induction-type tools have no sensitivity to resistivity anisotropy when 
die formation layers are perpendicular to the weU bore, but will be quite sensitive at 
relative dip angles above 60 degrees. See LOling et al.. Processing and modeling 2 MHz 
resistivity tools in dipping, laminated, anisotropic formaHons, Transactions of the SPWLA 
TmRTY-FlFiH ANNUAL LOGGING SYMPOSIUM. Tulsa. Oklahoma. June 19-22. 1994, paper 
QQ. 

Induction logging tools arc insensitive to resistivity anisotropy at low relative dip. 
The current density set up in the formation is limited to the horizontal direction. See 
Moran et al.. Effects of formation anisotropy on resistivity-logging measurements. 
GEOPHYSICS. V. 44. No. 7. July 1979. 1266-86; Anderson et al.. The Response of Induction 
Tools to Dipping Anisotropic Formations, Transactions of the SPWLA THmTY-SDOH 
ANNUAL Logging Symposium. Paris, ftance, June 26-29, 1995, paper D. At higher 
relative dips the measured resistivity will combine both vertical and horizontal 
components. However, the effect does not depend on tiie coil array spacing, so tiie 
components are not resolvable. See Anderson et al.. The Response of Induction Tools to 
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Dipping Anisotropic Formations, Transactions of the SPWLA Thirty-Sixth Annual 
Logging Symposium, Paris. France, June 26-29. 

Laterolog logging tools have some sensitivity to the vertical component of 
resistivity. See Moran et al.. Effects of formation anisotropy on resistivity-logging 
measurements, GEOPHYSICS, V. 44, No. 7, July 1979, 1266-86; Chemali et al„ The effect 
of shale anisotropy on focused resistivity devices. Transactions of the SPWLA TWENTY- 
Eighth Annual Logging Symposium, London, England, June 30-July 2, 1987, paper H. 
However, with the dual-Iaterolog, the deep laterolog was not veiy sensitive to anisotropy. 
The shallow laterolog, which returns current to nearby electrodes, is more sensitive to 
anisotropy {See Moran et al., 1979), but invasion will reduce the effect of anisotropy. 

It is desirable to obtain a reliable method and system for determining the resistivity 
properties and profiles of potential hydrocarbon-bearing zones in subsurface formations. 
Thus, there remains a need for improved methods and techniques to indicate formation 
resistivity anisotropy and evaluate the vertical component of resistivity. 

2. SUMMARY OF THE INVENTION 

Methods and systems are disclosed for combining induction tool and laterolog tool 
measurements to indicate that the formation resistivity is anisotropic and estimate the 
vertical component of resistivity even in the presence of invasion. 

One aspect of the invention provides a method for determining a characteristic of a 
subsurface earth fonnation surroundmg a borehole. The method comprises calculating at 
least one parameter representative of a property of the formation using log data obtained 
with a first well tool disposed within the borehole; modeling log data from the at least one 
calculated parameter; and comparing log data obtained with a second well tool disposed 
within the borehole against the modeled log data to determine the formation characteristic. 

Another aspect of the formation provides a method for determining a characteristic 
of a subsurface earth formation surrounding a borehole. The method cornprises disposing 
first and second well tools within the borehole to obtam measurement data; producing at 
least one parameter representative of a property of the formation from the measurement 
data obtained with the first well tool; producing theoretical data from the at least one 
produced panuneter; and comparing the theoretical data against the measurement data 
obtained with the second well tool to determine the formation characteristic. 

3 
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Another aspect of the invention provides a well logging system including fust and 
second well tools adapted to be moveable through a borehole and an apparatus adapted to 
be coupled to the well tools, the apparatus being adapted to respond to dau obtained with 
the well tools to determine a characteristic of a fonnation. The apparatus comprises means 
for calculating at least one parameter representative of a property of the formation using 
log data obtained with the first weU tool disposed within the borehole; means for modeling 
log data ftom the at least one calculated parameter; and means for comparing log data 
obtained with the second well tool disposed within the borehole against the modeled log 
data to determine the formation characteristic. 

Yet anotiier aspect of the invention provides a program storage device readable by 
a processor and encoding a program of instructions including instructions for performing 
operations in accord with the invention. 

3. BRIEF DESCRIPTION OF THE DRAWINGS 

Other aspects and advantages of the invention will become apparent upon reading 
tile foUowing detailed description and upon reference to tiie drawings in which: 

Figure 1 is a schematic diagram of tiie laterolog energization patterns obtained witij 
a conventional laterolog tool. 

Figure 2 is a schematic diagram of the current lines corresponding to the laterolog 
patterns of Figure 1. 

Figure 3 is a graph showing the response of the five laterolog anays of Figures 1 
and 2 in a homogeneous anisotropic formation. 

Figures 4-7 are graphical illustrations of the responses of tiie five laterolog arrays 
of Figures 1 and 2 in invaded, anisoti-opic formations. 

Figure 8 is a graph showing tiie comparison of a laterolog tool response in an 
anisotropic fonnation vs. its response in an isotiopic formation in accord witii tiie 
invention. 

Figure 9 illustrates several plots of a ratio used in tiie calculations of tfie invention. 
Figure 10 illustrates tiie ratio plots of Figure 9 witfi an applied scaling factor. 
Hgure 11 shows a crossplot for seventy ratios pertaining to cases of formation 
parametos in accord with the invention. 
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Figure 12 shows an average value of a parameter for the ratios of Figure 11 in 
accord with the invention. 

Figure 13 is a graph showing a computed formation vertical resistivity component 
in accord with the invention. 

Figure 14 is a graph showing computed induction tool and laterolog tool logs in 
accord with the invention. 

Figure 15 is a graph showing the inversion of the induction tool data of Figure 14 
using a laterolog model in accord with the invention. 

Figure 16 is a graph showing the vertical resistivity component derived from the 
data of Figures 14 and IS in accord with the invention. 

Figure 17 shows field logs obtained from an induction tool and a laterolog tool in 
accord with the invention. 

Figure 18 shows the vertical resistivity component derived from the data of Figure 
17 in accord with the invention. 

Figure 19 is a schematic diagram of a well logging system in accord with the 
invention. 

Figure 20 is a flow diagram of a method in accord with the invention. 
Figure 21 is a flow diagram of another method in accord with the invention. 

4. DETAILED DESCRIPTION OF SPEdHC EMBODIMENTS 

In the interest of clarity, not all features of actual implementation are described in 
this specification. It will be appreciated that although the development of any such actual 
implementation might be complex and time-consuming, it would nevertheless be a routine 
undertaking for diose of ordinary skill m the art having the benefit of this disclosure. 

Figure 1 shows the energization patterns, referred to as 'laterologs", obtainable 
from a laterolog-type tool 10 such as those disclosed in U.S. Pats. Nos. 5,754,050, 
5,852,363, 6,023,168, and Document EP-478409. The deepest array, mode 5, is 
considerably deeper than the LLS and should have considerable sensitivity to anisotropy. 
Figure 2 shows the current paths on several of the arrays of the tool 10. Note that the 
cunent flow lines turn parallel to the tool 10 axis and retum locally, thereby increasing the 
sensitivity to the vertical resistivity component. 

5 
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Figure 3 shows the response of the tool 10 arrays (1-5) to an anisotropic 
homogeneous formation. The ratio of the apparent resistivity Rapp over the formation 
horizontal resistivity component Rth is plotted against the ratio of the fonnation vertical 
resistivity component Rty over Rth. This separation shows that the shallow arrays are 
much more sensitive to an anisotropic formation than the deeper arrays. Although the 
homogeneous medium anisotropy response looks favorable, what really matters is the 
response in the presence of invasion. 

Figures 4 through 7 show four representative cases of formation resistivity 
measurements obtained with a laterolog tool in an anisotropic formation having invasion 
out to a radius of 60 inches. These figures show ratios of the formation horizontal 
resistivity component Rth over the formation invasion resistivity Rxo of 0.1, 0.2, 5, and 
10. For all the cases Rtv = 10 x Rth. The invaded zone is isotropic with resistivity Rxo. 
In Figures 4 and 5, Rxo > Rth. In Figures 6 and 7, Rxo < Rth. For all these cases, the 
invasion response confuses the anisotropy response. 

Figure 8 shows the anisotropic case of Figure 7 superimposed with the laterolog 
tool response to an isotropic virgin zone with formation resistivity Rt equal to Rth. As 
seen in Figure 8, the two curves representing Array 5 (RLA5) are almost parallel over the 
range of invasion radius out to 60 inches. The other shallower isotropic-vs.-anisotropic 
cases show separation at shallow invasion, but little or no separation at deeper invasion 
radii. The reasonable constancy of the difference between the isotropic and anisotropic 
Array 5 response suggests that the ratio of anisotropic to isotropic response or other ratios 
provide useful information. 

Designating Raa as the laterolog Array 5 log in an anisotropic formation, and Rai 
as the laterolog Array 5 log reading in an isotropic fonnation, the following ratio is 
derived: 

Rai 
(Raa -Rai)' 

Figure 9 shows a plot of the above ratio for Rxo =1 and Rth =10 for A values from 1.4 to 
10, where 
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It should be noted that the spread is ahnost proportional to X . This suggests that 
the ratio 

XRai 
Raa - Rai 



should be nearly constant. This ratio is shown in Figure 10. 



If the ratio 

A= (1) 
(Raa -Rai) 

is nearly constant for this case, then other values of Rxo and Rth should be looked at. 
Figure 1 1 shows a crossplot of 7 combinations of Rxo and Rth. There are 10 values of A 
ranging from 1.4 to 10. The spread in A is remarkably small. Figure 12 shows the average 
value of A from the data of Figure 1 1 . If an average value of A is taken to be a constant (or 
fitted function of the radius of invasion n), then the unknown A can be solved for 

Rai 

One unknown remains in equation (2) - Rai. 

As mentioned above, induction tools have litfle or no sensitivity to Rtv in a vertical 
well. In such a case» the mduction tool log data is inverted for Rxo, Rth, and ri using a 
lookup table, computer algorithm, model, or any other inversion means as known in the 
art. Also known in the art is the use of modeling techniques to simulate specific tool 
responses from given data. By applying a one-dimensional (ID) model of the laterolog 
tool response to the formation parameters obtained from the induction tool, Rxo, Rth, and 
ri, the laterolog log data in an isotropic foimation (Rai) can be inferred. 

Then equation (2) can be solved for X , and then 

Rtv=A^Rth. (3) 
Figure 13 shows such a solution using a fitted value of A as a fiinction of ri. The values 
for A may be precomputed. The value for A ID model for the laterolog tool was derived 
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using a simple average 7-function. which is only a gross estimate for an example. The 
estimate of Rtv is within a factor of 2 of the actual value of Rtv in the uninvaded zone. 
Although a factor of 2 may be considered inaccurate for estimates of petrophysical 
quantities, such accuracy is ample for determining that anisotropy exists in a vertical well. 
The accuracy may be improved by refining A. Both induction and laterolog tools have 
effects in 2D formations (layers and invasion in the layers) that are not accurately modeled 
with ID forward models. To test the usefulness of the ID models used in the recreation of 
Rai, it is advisable to try the algorithm in a 2D-formation model. 

Figure 14 shows computed induction and laterolog tool logs in a simple 2D 
formation with 3 invaded, anisotropic layers of varying thickness. These layers have 
identical values of Rtv. Rth. Rxo, and n. The thickness varies ftom 50 ft. to 10 ft. The 
induction tool log data reflects the receiver to transmitter spacings (10 in.- 90 in.) for the 
particular tool. 

Figure 15 shows a modeled laterolog tool log derived from the induction tool log 
data of Figure 14. The modeled log data represents the theoretical log values obtainable 
with the laterolog tool disposed within the isotropic formation. The modeled laterolog 
log data are tiien compared against the actual or apparent laterolog log data to compute 
Rtv as described above. Figure 16 shows the resulting Rtv log computed from the data of 
Figures 14 and 15. 

The Rtv log of Figure 16 shows anomalies in the non-invaded isotropic beds 
between the invaded anisotropic beds. As seen in Figure 14, the laterolog tool logs 
display shoulder effect in these beds and on eitiier side of the first and third beds. Any 
source of error that leads to an increased separation between the laterolog tool logs 
predicted from the induction tool parameters and apparent laterolog tool logs wiU be 
interpreted as anisotropy. OveraU, however, intelligent interpretation should remove 
most ambiguity. 

A field test was conducted comparing a laterolog tool against an induction tool in a 
formation with expected resistivity anisotropy. Figure 17 shows the induction tool and 
laterolog tool logs in a pay sand. Figure 18 shows the indicated vertical resistivity Rtv 
along with the induction tool Rth, Rxo. and ri inversion in a zone that exhibits the 
anisotropy signature. It also exhibits shoulder effect just as in the computed example of 
Figure 16. 
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4.1 Program Storage Device 

It will be apparent to those of ordinary skill having the benefit of this disclosure 
that the present invention may be implemented by progranuning one or more suitable 
general-purpose computers having appropriate hardware. The progranraiing may be 
accomplished through the use of one or more program storage devices readable by the 
computer processor and encoding one or more programs of instructions executable by the 
computer for perfonning the operations described above. The program storage device may 
take die form of, e.g., one or more floppy disks; a CD ROM or other optical disk; a 
magnetic tape; a read-only memory chip (ROM); and other fonns of the kind well-known 
in the art or subsequently developed. The program of mstructions may be "object code," 
i.e.. in binary form that is executable more-or-less directly by the computer; in "source 
code" that requires compilation or interpretation before execution; or in some intermediate 
form such as partially compiled code. The precise forms of the program storage device 
and of the encoding of instmcdons are inunaterial here. 

42 Logging System 

The present invention can be in^)lemented in a logging system 200 with a first well 
tool 210 and a second well tool 215 adapted to be moveable through a borehole and an 
apparatus 220 coupled to the tools as shown in Figure 19. Certain conventional details are 
omitted in Figure 19 for clarity of illustration. The apparatus 220 comprises a surface 
computer (such as the general-purpose computer and program storage device described 
above) coupled to the well tools 210, 215 by a wirelme cable 225 or linked to die tool 
string's surface instrumentation as known in the art. 

The first well tool 210 may be any laterolog-type tool that returns current locally, 
i.e., to the sonde body. The second well tool 215 may be any propagation or induction tool 
that provides multiple curves, preferably providing at least three independent 
measurements to solve for Rxo, Rt, and n. In addition to the U.S. Patents mentioned 
above, other useable tools include the instruments described and illustrated in U.S. Pats. 
Nos. 5,339,037 and 4,968,940, bodi assigned to die present assignee. 

The formation parameters can be determined in real-time by sending die 
measurement data to the surface as they are acquired, or it can be determined from a 
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recorded-mode by recording the data on a suitable recordable medium. As known in the 
art, the measured data are transmitted from the tools 210. 215 to the surface computer 220 
by electronics (not shown) housed in the tools. The data may be sent to the surface 
computer 220 along the wireline cable 225 or by mud-pulse telemetry if implemented in a 
logging-while-drilling configuration. It will be understood that alternative means can be 
employed for communicating the acquired data to the surface as the precise forai of 
conomunication is immaterial to the implementation of the disclosed techniques. 

Once received by the surface computer 220, the data can be recorded, processed, or 
computed as desired by the user to generate the formation resistivity profile. The profile 
can then be recorded on a suitable output record medium. Alternatively, some or all of the 
processing can be perfonned downhole and the data can be recorded uphole. downhole, or 
both. It will be understood that the present invention can be implemented in conjunction 
with any suitable technique for keeping track of the tool 210, 215 depth within a borehole. 

Figure 20 illustrates a flow diagram of a method for determining a characteristic of 
a subsurface formation surrounding a borehole in accord with the invention. The metiiod 
comprises calculating at least one parameter representative of a property of the formation 
using log data obtained witii a first weU tool disposed witiiin the borehole 300; modeling 
log data from tiie at least one calculated parameter 305; and comparing log data obtained 
wifli a second well tool disposed within the borehole against the modeled log data to 
determine the formation characteristic 310. 

Figure 21 iUustrates a flow diagram of a metiiod for determining a characteristic of 
a subsurface foimation sunounding a borehole in accord with the invention. The metiiod 
comprises disposing first and second weU tools within tfie borehole to obtain measurement 
data 400; producing at least one parameter representative of a property of tiie formation 
from tile measurement data obtained witii tfie first weU tool 405; producing theoretical data 
from tiie at least one produced parameter 410; and comparing the produced tiieoretical data 
against tiie measurement data obtained witii the second well tool to determine flie 
formation characteristic 415. 

WhUe tiie metiiods and apparatus of tiiis invention have been described as specific 
embodiments, it will be apparent to tiiose skiUed in tiie art tiiat variations may be applied 
to tiie structures and in the steps or in tiie sequence of steps of tiie metiiods described 
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herein without departing from the concept and scope of the invention. All such similar 
variations apparent to those skilled in the art are deemed to be within this 
concept and scope of the invention as defined by the appended claims. 
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CLAIMS 

1. A method for determining a characteristic of a subsurface earth formation 
surrounding a borehole, comprising: 

(a) calculating at least one parameter representative of a propeny of the 
formation using log data obtained with a first well tool disposed within 
the borehole; 

(b) modeling log data from the at least one calculated parameter; and 

(c) comparing log data obtained with a second well tool disposed within 
the borehole against the modeled log data to deteraiine the formation 
characteristic. 

2. The method of claim 1, wherein step (b) comprises modeling log data 
representative of log data theoretically obtainable with the second well tool 
disposed within the borehole. 

3. The method of claim 1 , wherein step (c) comprises calculating a ratio between 
the log data obtained with the second well tool and the modeled log data. 

4. The method of claim 1, wherein step (c) comprises calculating the difference 
between the log data obtained with the second well tool and the modeled log 
data. 

5. The method of claim 4, wherein step (c) comprises multiplying the calculated 
difference by a predetermined factor. 

6. The method of claim 1, wherein the first well tool is an induction-type tool and 
the second well tool is a laterolog-type tool. 

7. The method of claim 1, wherein the determined fomiation characteristic is a 
resistivity value or resistivity profile and the method is performed during or 
after drilling of the borehole. 
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8. A method for determining a characteristic of a subsurface earth formation 
surrounding a borehole, comprising: 

(a) disposing first and second well tools within the borehole to obtain 
measurement data; 

(b) producing at least one parameter representative of a property of the 
formation from the measurement data obtained with the first well tool; 

(c) producing theoretical data from the at least one produced parameter of 
step (b); and 

(d) comparing the theoretical data of step (c) against the measurement data 
obtained with the second well tool to determine the formation 
characteristic. 

9. The method of claim 8, wherein step (c) comprises producing theoretical data 
representative of measurement data obtainable with the second well tool 
disposed within the borehole. 

' 10. The method of claim 8, wherein step (d) comprises calculating a ratio between 
the measurement data obtained with the second well tool and the produced 
theoretical data. 

11. The method of claim 8, wherein step (d) comprises calculating the difference 
between the measurement data obtained with the second well tool and the 
produced theoretical data. 

12. The method of claim 1 1. wherein step (d) comprises multiplying the calculated 
difference by a predetermined factor. 

13. The method of claun 8, wherein the first well tool is an induction-type tool and 
the second well tool is a laterolog-type tool. 

14. The method of claim 8, wherein the determined formation characteristic is a 
resistivity value or resistivity profile and the method is performed during or 
after drilling of the borehole. 

13 



BNSOOCIO: <GB_2afle9lSA_l_> 



A well logging system including first and second well tools adapted to be 
moveable through a borehole and an apparatus adapted to be coupled to the 
well tools, the apparatus being adapted to respond to data obtained with the 
well tools to detennine a characteristic of a formation, the apparatus 
comprising: 

means for calculating at least one parameter representative of a property of 
the formation using log data obtained with the first well tool disposed 
within the borehole; 

means for modeling log data from the at least one calculated parameter; 
and 

means for comparing log data obtained with the second well tool disposed 
within the borehole against the mo(teled log data to detennine the 
formation characteristic. 



The system of claim 15, the modeling means further comprising means for 
modeling log data that is representative of log data theoretically obtainable 
with the second well tool disposed within the borehole. 

The system of claun 15, the conq)arison means further comprising means for 
calculating a ratio between the log data obtained with the second well tool and 
the modeled log data. 

The system of claim 15, the comparison means further comprising means for 
calculating the difference between the log data obtained with the second well 
tool and the modeled log data. 

The system of claim 15, wherein the first weU tool is an induction-type tool 
and the second well tool is a laterolog-type tool. 

The system of claim 15, wherein the determined formation characteristic is a 
resistivity value or resistivity profile. 
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